The critical Ising lines of the d = 2 Ashkin- Teller model 
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The universal critical point ratio Q is exploited to determine positions of the critical Ising 
transition lines on the phase diagram of the Ashkin- Teller (AT) model on the square lattice. A 
leading-order expansion of the ratio Q in the presence of a non- vanishing thermal field is found from 
finite-size scaling and the corresponding expression is fitted to the accurate perturbative transfer- 
matrix data calculations for the L x L square clusters with L < 9. 



The AT model has first been proposed as a model of a 
four component alloy (j]]. It has attracted a lot of theo- 
retical interest for years because it is a simple and non- 
trivial generalization of the Ising and four-state Potts 
models. Fan g] has shown that the hamiltonian of the 
AT model can also be written with two Ising variables 
(S = ±l,cr = ±1) located at each site of the lattice, 
which in the presence of a magnetic field has the form: 



T~t = — ^ (JiSiSj + J2<Ji&j 

<i,j> 

N 

+ J4Si<TiSj(Jj + J ) — h'y^SjQ-j 



(1) 



Herein we consider only the nearest neighbour pair inter- 
actions on the simple square lattice consisting of N = L 2 
sites with periodic boundary conditions and we assume 
that Ji = J2 (isotropic case). 

Wagner H has shown that the AT model is equivalent 
to the alternated eight vertex model, which has not been 
solved exactly. Only one critical line in the phase diagram 
of the isotropic AT model is known exactly thanks to the 
duality relation found by Fan For this reason many 
approximate approaches have been applied for construct- 
ing the complete phase diagram: the mean field theory 
(MFA) mean-field renormalisation group (MFRG) 

M, renormalisation group (RG) 
renormalisation group (MCRG) 
paper to establish an accurate location of the remaining 
critical lines. 

In our approach we exploit finite-size scaling for the 
ratio of the square of the second moment to the fourth 
moment of the order parameter M: 



and Monte Carlo 
It is the aim of this 



<M 2 >| 
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where < ... > means thermal average and the index L 
indicates the linear size of the system (L x L). In the 
limit L — > 00 this ratio becomes universal in the critical 
point |hJ and is denoted Q hereafter. Three not ex- 
actly known critical lines of the isotropic AT model are 
believed to belong to the Ising universality class [p|,pd|. 
Here it is assumed that these lines correspond to the 
Ising-like continuous transitions with the order param- 
eter M = J2iLi Si<Ji. A scaling formula for Ql can be 
derived starting from the finite-size scaling relation for 
the singular part of the free energy for the square Ising 
model 



{gug^L- 1 ) = A(g t L) InL + B(g t L, g h V">) (3) 

where A and B are unknown amplitudes, gt, gt are non- 
linear scaling fields and is the magnetic critical ex- 
ponent. The nonlinear scaling fields gt and g^ can be 
expanded in terms of the corresponding linear thermal 
and magnetic scaling fields t and h. 

Taking into account the relations between the magneti- 
zation moments in Eq. (0) and the corresponding deriva- 
tives of the free energy we have calculated the scaling 
expansion for Qi(i, h = 0) to the leading order in t and 
up to L 3 ~ 4yh : 



M = QUO) 



dQ L (t) 



dt 



t 



(4) 



The zcroth order term <5l(0) was evaluated previously 
[fl2| and the first order term is of the form 



1 



dQ L (t) 



dt 



= aiL + a 2 + a 3 L 3 - 2yh + (a 4 + a 5 \D.L)L 2 ~ 2yh + a 6 L 5 - 4yh + 

t=o 

+ (a 7 + a 8 XnL)^- 21 ^ + (a g + a w In L)L 4 - iyh + a n L- 2yh + 
-(au + a 13 In L)L 7 - &Vh + (a u + a w In L + a 16 In 2 L)L 3 - iyh + ... , 



(5) 



where on (i — 1, ...16) are unknown amplitudes. In our 
work we consider only the first three terms in the expan- 
sion (§), but for some future Monte Carlo applications 
the higher order terms in 1/L might be important. 

We have calculated the QL(t) ratio exploiting the 
transfer matrix technique which for the Ising model was 
explained in . Our system consists of L columns con- 
taining L sites. Spins from the jth column are denoted 
by Ej = (Sji,<jji, Sj 2 ,CTj2, —,Sj L ,o-j L ) so that 



Z= ex P(~/3#(Ei, Et)) = TtT 1 



(6) 



Si ,£2, 



where T is a 4 L x 4 L transfer matrix. This can be split 
into the product T = T^T W of a diagonal matrix TV and 
a non-diagonal matrix T/ t containing the intra- and the 
inter-column interactions respectively. They are defined 
as follows 



( N \ 

T„(Efc, E;) = <5g fc g i exp I y^^SkjSkj+i + K 2 a kyi (Tk,i+i + K4Sk,i&k,iSk,i+i<rk,i+i + HSk,i<Tk,i) J 



(7) 



T^(E fc , S;) = exp y^{K 2 Sk,iSi. l + K 2 <rk,i<ri,i + K 4 Sk^ak,iSi ti ai t 
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where , = -^ T , K t = J t (3 (i = 1,2,4) and H = f3h. 
The latter matrix can be expressed as a product of sparse 
matrices which facilitates the numerical calculations. 

The averages in Eq. (|^) can be expressed in terms of the 
corresponding coefficients Zk ]l2| in the expansion of the 

field dependent partition function Z(h) = X^o-^fclZ- 
The coefficients Zk can then be calculated from Eq. (g) 
by multiplying the base vectors by matrices T„ and T/j 
in such a manner that the terms in the same power of h 
are kept separately jTJj . 

At first we calculate the amplitudes a>i (i < 5) from Eqs 
(0) and (^) with known values Ql(0). In the limit K 2 — 
0, i.e. the Ising model in Sa, K± c — K c — |ln(l + y/2) 
and in this case we have only one coupling constant (-K4). 
Thus we can write the reduced temperature in the form: 



K, 



4c 



K 4 



K 4 



(9) 



Selecting different values of the scaling field t we can 
solve the set of linear algebraic equations for a^. For the 
ferromagnetic coupling K4 we consider the system sizes 
L = 2, 3, 9 whereas for the anti-ferromagnetic one only 
the even values L — 2, 4, 6, 8 are considered, so that we 
can evaluate the coefficients on up to i = 5 or i = 3, 
respectively. 

Having fixed K 2 ^ and knowing the en (i < 3) and 
Ql(0), we have calculated Ql{K 2 ,K/C) for a number of 



couplings K 4 . This enables a determination of the corre- 
sponding t values from Eqs (|4|) and (5). Then knowing 
t we can easily obtain K 4c from Eq. (3) and K 2c from a 
similar equation, but written for K 2 . The estimates K 4c 
and K 2c are very stable if we find t e< 1CP 7 , 10~ 4 >. 

The exactly known critical curve with continuously 
varying critical exponents fll|] is terminated in the 4- 
state Potts point where it bifurcates. In the vicinity of 
this point the convergence of our results is diminished 
and the estimates of K 4c become size dependent. This 
size dependence is illustrated in Fig 1. Due to the limited 
number of system sizes available in our calculations we do 
not try to include any corrections to scaling and we sim- 
ply extrapolate our data. The corresponding estimates 
are shown on the ordinate axis in Fig 1. Such a strong 
size dependence does not occur for the anti-ferromagnetic 
couplings, since there is no Potts point in this case. 

Our final results represented by open circles connected 
by thin continuous lines are shown in Fig. || and they are 
compared with other results and predictions. The numer- 
ical uncertainties do not exceed the size of the symbol. 
The curve plotted by the bold line represents the part of 
the phase diagram found exactly by Baxter jll| . It sep- 
arates the Baxter phase B from the paramagnetic phase 
P. The ferromagnetic and anti-ferromagnetic phases with 
non-vanishing order parameter M are denoted by the la- 
bels F and AF, respectively. 

In the ferromagnetic region K 4 > we have only cal- 
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culated the curve joining the 4-state Potts point to the 
pure Ising point K c at K2 = 0. The second branch fol- 
lows from the corresponding duality relation |5||i~i|]. In 
the boundary between AF and P phases with the dotted 
lines we plot the approximate curve as given by Baxter 
fllf and in the ferromagnetic region we also include the 
MCRG results marked by filled circles. 

As can be seen (Fig. ||) our results are in good 
agreement with the MCRG || approach, but are quite 
different from Baxter's predictions |Tl| in the anti- 
ferromagnetic region. For the boundary between AF and 
P phases, our results coincide with those obtained by 
Mazzeo et al. |l3| |. These authors actually investigated 
the six vertex model with the transfer matrix technique in 
combination with conformal invariance arguments; their 
results can be mapped onto the results for the P-phase 
boundaries and they are shown in Fig. 2. 

As to our accuracy: near the ferromagnetic Ising point 
it is of about 2 * 1CP 6 and in the neighbourhood of the 
Potts point it decreases down to about 3 * 10~ 2 . The 
accuracy in the anti-ferromagnetic region is even better: 
near the Ising point it reaches 5 * 10~ 8 and for the high- 
est point at the phase diagram in Fig. || it decreases to 
3* 1(T 3 . 
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FIG. 1. The L-dependence of the critical values of the parameter K^ c . The points on the vertical axis are the extrapolated 
values. 
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FIG. 2. Phase diagram of the AT model in two dimensions. The solid bold line represents the exactly known critical line, 
which terminates at the 4-state Potts point. Empty circles with continuous lines describe our results. The solid circles display 
MCRG results , dotted line is drawn after Baxter and diamonds are the transfer matrix results combined with conformal 
invariance. 
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